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ABSTRACT 

An accurate theoretical description of structure formation at least in the mildly non-linear 
regime is essential for comparison with data from next generation galaxy surveys. In a recent 
approach one follows the time evoluti on of correla tors directly and finds a hierarchy of evolu- 
tion equations with increasing order ( Pietronil2008 1. So far, in this so called time renormalisa- 
tion group method the trispectrum was neglected in order to obtain a closed set of equations. 
In this work we study the influence of the trispectrum on the evolution of the power spectrum. 
In order to keep the numerical cost at a manageable level we use the tree-level trispectrum 
from Eulerian perturbation theory. In comparison to numerical simulations we find improve- 
ment in the mildly non-linear regime up to k 0.25 hMpc . Beyond k 0.25 hMpc the 
perturbative description of the trispectrum fails and the method performs worse than without 
the trispectrum included. Our results reinforce the conceptual advantage of the time renormal- 
isation group method with respect to perturbation theory. 



1 INTRODUCTION 

In the contemporary picture of our Universe structures evolve from nearly Gaussian distributed small perturbations in the homogeneous 
density field. Sound waves formed in the coupled photon-baryon fluid before recombination left oscillatory features in the matter power 
spectrum, so called baryonic acoustic oscillations (BAO). Detections of this efFect have become a valuable too l to constrain cosmo l ogical 



parameters since ampl itude and position of the o scilla tions depend on the expansion history of the Universe ( Eisenstein et al. 2005; 



Hiitsi 



2006; Eisenstein et al 112007 ; Padmanabhan et al. 2007; Blake et al. 2007). Current and upcoming galaxy surveys - such as BQSS3,WFIRST 



HETDEX Mill et al 



2008) and WFMOS iGlazebrook et al. 2005) - will measu re the power spectrum of the matter d istribution to 1% 



a 

accuracy in the region of the baryonic oscillations 0.05 ZiMpc 1 < k < 0.25 AMpc 1 l Eisenstein et alfl998 | ; Seo & Eisensteinll2003 ). 

Since mode-coupling effects can significantly influence the position of the first peak at low redshifts JCrocce & Scoccimarrol2008 ), it is 



necessary to find a robust theoretical description of structure formation in the linear and mildly non-linear regime. While standard perturbation 
theor y l Bernardeau et al 1l2002h is a powerful tool for comparison with ob servations from galaxy surv eys on large scales Jjeong & KomatsJ 
2006t 20091) . it breaks down at the scales of baryonic acoustic oscillations 1 Jain & Bertschingerlll994 ). 



The i ncrease of computational power and efficiency of algorithms made iV-body simulations the most established approach f or structure 



formation JSpringelkooltlHuff et alfeoollHeitmann et alkoojjTakahashi et alj2008l ; lEvrard et al.l2008l;lHeitmann et al.l20icl) . However 
to extract statistical information from numerical simulations, either large sets of initial conditions or large volumes are needed and it is 
difficult to control measurement uncertainties in the mildly non-linear regime ( lAngulo et al.l 120081) . In a recent work a simple physically 
motivated picture was used to reduce the sample variance to speed up the scanning for cosmological parameters with /V-body simulations 
1 Tassev & Zaldarriagalboi 1 ), but it remains numerically expensive and we favour a fast semi-analytical tool to predict non-linear structure 

formation. 

While the ha lo model approach 1 Peacock & Doddsll996l ; Smith et al .120031) was found to be uncapable of reaching the required accuracy 
(Huff et al. 2007), different attempts to include corrections of specific types to all orders at the same time have been presented ove r the 
last years. Field theoret ical techniques motivated by t he renormalisation group jMatarrese & PietronilEo07l I2OO8I lAnselmi et al]|201 lb and 
resummation methods 1 Crocce & Scoccimarroll2006allb ) improved the results for power spectra in the mildly non-linear regime sig nificantly 
down to z = in comparison to /V-body simulations. In the following also higher order statistics were studied in these frameworks ( Valageas 



1 http://www.sdss3.org/surveys/boss.php 

2 http://wrirst.gsfc.nasa.gov 
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2008; Bernardeau et al .|2008| ; IGuo & Jinj2009l) . However, these approaches are formulated for an Einstein-de Sitter cosmology and are later 
generalised to other cosmologies by s ubstituting the respective growth function. The accuracy of this approximation is hard to quantify at 
higher orders dBernardeau et al 1120021) . 

Correlating the structure formation equations with fields to different orders leads to an in finite hierarchy of evolution equations for 
correlators with increasing order. This hierarchy - si milar to the we ll known BBGKY hierarchy dPeeblesll 19801) - was truncated at the level 
of the trispectrum to obtain a closed set of equations ( Pietronill2008 ). This so called time renormalisation group method (TRG) can easily be 
generalised to a large set of di fferent cosmologies including models with scale dependent growth functions. For example, this i s important in 
the case of massive neutrinos dLesgourgues & Pa stor 2006) which was also studied within this framework dLesgourgues et al ]|2009l) . 

The objective of this work is to study the effect of the trispectrum on the TRG approach. Including the entire time evolution of the 
trispectrum would lead to an immoderate numerical effort. We therefore include the tree-level trispectrum from perturbation theory into 
the evolution equation of the bispectrum and study its effects on the power spectrum in the BAO regime. We work in a standard ACDM 
cosmology close to the best-fitting cosmology (£1 M = 0.25, fl b /i 2 = 0.0224, h = 0.72, n = 0.97 and <r 8 = 0.8) and compare to power spectra 
obtained from Af-body simulations of the same cosmology dCarlson et al . 2009). 

The paper is organised as follows: In Section [2] we review the time renormalisation approach introduced by 1 



Pictroni 



(120081) . How the 

perturbation theory trispectrum can be included into this method will be subject of Section[3] In Section|4]we include the trispectrum into the 
formal analytic solution of the system and discuss the additional corrections in a diagrammatic representation and in Section[5]our numerical 
results and their comparisons to A'-body simulations are presented. The results are summarised and discussed in Section|6] 



2 TIME RENORMALISATION 



In thi s section we write the structure formation equations in a compact matrix form and review the time renormalisation approach dPietroni 



2008) as a starting point for our further calculations. We will concentrate on spatially flat cosmologies with a dark matter component and a 



non-clustering dark energy fluid or ACDM. How ever, this metho d can be easily extended to more exotic cosmologies including cosmologies 
with a scale-dependent growth function D + (k, a) dPietronill2008h . 



2.1 Non-linear Structure Formation 

The time evolution of the dark matter density contrast S(x, t), the peculiar velocity v(x, t) and the fluctuation of the gravitational potential 
0(x, t) is governed by the continuity-, Euler- and Poisson equations: 

<9 r <5(x) + V-[(l+<5(x)) v(x)] = (1) 
d T v(x) + <H v(x) + (v(x) ■ V) v(x) = -V <D(x) (2) 

V 2 (D(x) = ^<H 2 Sl m 6(x). (3) 

Here, the Hubble function "H is the logarithmic derivative of the scale factor with respect to conformal time: "H = dlna/dr. For notational 
clarity the time dependences will be omitted for most of the quantities throughout the whole work. In the case of a one-component dark 
matter fluid we have fl m = 1, while in presence of an additional dark energy fluid we only have to replace Sl m by 

o.-(i + f«-'f. 

The quantities p° is _ and pjj, represent the background densities of the dark energy and dark matter fluid components respectively, w denotes the 
equation of state parameter of the dark energy fluid, p dc = wpL. Models with a cosmological constant, as for example in ACDM, can then 
simply be described by w = -1. We assume a irrotational peculiar velocity field, V x v(x) = 0. Due to Helmholtz's theorem the velocity field 
can th en be described by its d ivergence only, 9(x) = V • v(x). For the two remaining scalar fields 9(x) and 6(x) the evolution in Fourier space 



reads (Bernardeau et al. 2002): 



)ed by 1 
et alb 



d T S(k) + 9(k) = -<5 D (k - q - p) a(q, p) <9(q) <5(p) 
d T 6(k) + 'He{k)+l'H 1 il m S(k) = -<5 D (k-q-p)£(q,p)0(q)0(p). (5) 

Here and in the following, repeated momentum arguments in products (for instance p, q in eqn. lO) imply an integration over the respective 
momenta f d 3 p j d 3 q... . The left hand side of eqn. {5]l represents the linear evolution of the fields while the right hand side expresses the 
non-linear mode-coupling, which is determined by the two model independent mode-coupling functions 

(p + q) q 

a(q,p) = , 

T 

2q 2 p 2 
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To write the equations in a more compact form, we introduce the logarithmic time variable 77, which explicitly contains the linear structure 
growth with respect to the scale factor: 



The scale factor a m may be chosen at a time when the systen could still be well aproximated to be Gaussian and in the linear regime. Now a 
doublet field can be introduced for the density contrast c5(k) and the divergence of the velocity dispersion 9(k): 



v ,(k) \ _„ I <5(k) 
I 6 I -fl(k)/(W/ + ) 



(8) 



n «b-\ lfn4 as on I ' (10) 



Here, /+ is the logarithmic derivative of the growth function D + with respect to the scale factor a, f + = dlnD+/dlna. The factor e - '' now 
compensates for the linear evolution of the fields <ft(k). In other words, solving the linearised structure formation equations for these fields 
would lead to no time dependence in the doublet field at all. Therefore, any evolution away from the initial field douplet is explicitly due to 
non-linear effects. The set of equations ((5) can now be expressed in a very compact form: 

d,, <^ n (k) = Q. ab <p b (k) + & y abc (k, -q, -p) <p b (q) tp c (p) . (9) 
Here, the linear evolution is governed by the matrix 

1 -1 

-3fi m /(2/ + 2 ) 3fi m /(2/ 2 ) 

while the non-linear mode-coupling is moderated by the vertex funtions y abc (k, q, p). The only non- vanishing vertex functions are 

a (q, p) 

7121 (k,q,p) = y 112 (k,p,q) = (5 D (k + q + p) — - — , 

y 2 22 (k,q,p) = c5 D (k + q + p)/3(q,p) . (11) 

Since the vertex functions only appear in integrals, it is notationally convenient to introduce vertex functions, for which the ^-function is 
already integrated out: 

y*** s J dVy*(k,q,p). (12) 
These quantities turn out to dependent only on the absolute values of the wave vectors k, q and |k + q|. 



2.2 Hierarchy of Correlators 

While in standard perturbation theory one aims to solve the evolution equations of the fields themselves, in the time renormalisation approach 
one formulates evolution equations directly for the final quantities of interest - the correlators of the fields. In order to do this, one can use 
eqn. {9]l to write down a hierarchy of evolution equations for correlators of any order: 

dq(<p a <p b ) = -Cl ac {ip c (pi,) - Cl bc ((p a <p c ) + el [facditPctPdVb) + fbcd((Pa<Pc<fid)] 

d n (tp a <p b <Pc) = -QadiVdVbVc) + c.p.fa, b, c] + e'' [yadeiVdVeVbVc) + c.p.fa, b, cj] 

d,,(<Pa(PbtPc<Pd) = ••• • (13) 

Here and in the following, we abbreviate the doublet field index a t and the wave vector ki to a single number index, i.e. tpi = ^> Hl (ki). As 
a natural property of this hierarchy, for the evolution of a correlator of order n the knowledge of correlators of the next higher order n + 1 
is needed. Therefore, one is obliged to truncate this hierarchy at a certain point in order to obtain a closed set of equations. Splitting up the 
four-point correlator into its connected and unconnected part yields by Wick's theorem 

(<Pl<P2<P3<P4> = (<Pl<P2</>3<P4>c + (VlVzXVdV*) + (.fm)(.V2f4) + (Vl^iVzVi) ■ ( 14 ) 

For instance, one can close the system by neglecting the connected part of the four-point correlator, as it was done by Pietroni 1 2008h . Due 
to this approximation, one is left with the first two equations of the hierarchy and the simplified system is then fully described by its power 
spectra PK and its bispectra B ''* 2 '* 3 

W 2 > = (5 D (k, +k 2 )^ 

(<pm<P3) = S D (k, + k 2 + k 3 ) B k J* 2 ' h . (15) 

Due to isotropy the bispectrum will only depend on the absolute values of the wave vectors. Integrating eqn. d 1 3 fc over one wave vector and 
using eqn. U4b with {<pi<P2<P3'P4)c = one finds a closed system of equations in which the four-point function is represented in terms of power 
spectra P\: 



fl P* 

U ri r ab 
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Here and in the following p will denote the vector p = -(k + q). How this system can be solved numerically will be subject of the next 
subsection. 



2.3 Solving the Closed System 

Formal and numerical solutions to the closed system in eqn. ( I16t have already been presented ( Pietronil2008 ). We follow the same numerical 



7* 

acd,bef 



path to investigate the solution's sensitivity to the perturbation theory trispectrum. Since we are mainly interested in the evolution of the 
power spectrum itself, one circumvents the necessity of tracking the total bispectrum by introducing auxiliary integrals 

The introduction of these integrals encapsulates the one loop character of the power spectrum evolution equation (first equation in l|16t), 
which simplifies to the tree level equation 

V* = - Ofc^i + e"^ [C dM . d + lU cd ] . (18) 

Differentiating the integrals l k cdbej - with respect to time 77 and using the evolution equation of the bispectrum from eqn. l l 1 6t gives the 
following time evolution of these integrals: 

^1^acd,bef = ~QbgI a cd,gef ~ ^S^acdjjgf ~ ^fg^acdjxg + ^ acd, bef 1 (^) 

The £-space loop integral which was originally in the evolution of the power spectrum now appears in the mode-coupling integrals k k acdbe p 
which initially drive the system away from Gaussianity: 



acdMef 



^ J * 3 4 W (y$X p l pP M + yT p lA + y'KA) + (* ~ p)} ■ ( 20 > 



The calculation of the integrals A^ cdb . is so far the only time consuming task in numerically solving the system. 
2.4 Initial Conditions and Symmetries 

We propagate the system of equations <18M19t forward in time starting from an initial time at which the dynamics of the fields could still 
be well approximated by the linearised evolution equations. As initial redshift we choose z in =100 and start with the linear power spectrum 
as initial contitions: 

P k ab (r, = 0) = u a u b P£(77 = 0) , (m.m) = (1, 1) . (21) 

Furthermore, we will assume Gaussian initial conditions. This implies a vanishing inititial bispectrum B*£? and therefore also vanishing 
initial integrals Mf : 

iZiMfb = °> = KZ P( V = 0) = . (22) 

Following the next argments, the full system of 64 integrals l k cdbe j and 3 power spectra P k ab can be reduced to 14 independent components. 
Since the only non-vanishing vertex contributions appear for the index triples (acd) € {(112), (121), (222)) only the integrand I k cdbe j, with 
these triples will evolve away from zero, which can be seen from equations l|19l>-J20t. The remaining 24 components can be further reduced 
by symmetry arguments. Using the following symmetries in the vertex functions and the bispectrum 

7^7=7^", Bt**=B%', (23) 

we can find the following symmetry for the integrals I k cdbe f- 

I acd, be f = Kdcbfe ■ (24) 

Due to this symmetry, only 14 independent integrals remain to be followed. The independent integrals are identified by the direct product of 
(acd) = (112) and (bef) = (M1),(M2), (£21) and (622), (b = 1,2), and the direct product of (acd) = (222) and (bef) = (M1),(M2) and 
(b22), (b = 1, 2). Including the 3 independent power spectra P* p Pf, and PL implies a system of 17 components in total. 



3 TRISPECTRUM 

The connected part of the four-point correlator is called the trispectrum. To investigate the method's sensitivity with respect to the trispectrum, 
we include the perturbation theory trispectrum to third order in the linear power spectrum P^. Including the full time evolution of the non- 
perturbative trispectrum would increase the numerical effort disproportionally. This would be given by the third equation in the hierarchy 

© 2012 RAS, MNRAS OOO.mfTol 
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of time evolution equations for correlators in eqn. (113) . However, the time evolution of the tree-level perturbation theory trispectrum is 
exclusively given in terms of the growth function D + (a)l D + (a m ) = e 17 . Therefore, using this approximation for the bispectrum it is sufficient 
to compute trispectrum corrections at the initial time a m . These corrections can then be included into the routine without increasing the 
computational cost significantly. 



3.1 Perturbation Theory 

In the following we adapt standard perturbation theory dBernardeau et al .120021) to the compact matrix formulation of structure formation. We 
expand the fields <p a (k, rf) in «th-order perturbative contributions y>„ (k, rf), which can be written in terms of the perturbation theory kernels 
F< n) (ki, k„) and the initial linear fields ^"(k) = ^ fl L (k, rj = 0): 

<p a (Kv) = J^'rf'M (25) 

jp£°(k) = fdV- fd 3 9 „5 D (k-qi... B )FW(q 1) ... ) q n )^ 1 >(q 1 )...^ I) (q n ) (26) 

with qi „ = qi + ... + q„. For notational simplicity we combined the two standard kernels of each order into a vector, since we will also need 
trispectrum correlations to the velocity field components: 

7 <„)„. ... / ^ n) (ki,...,k„) 



/,k k i \ <r ^..:^i- (27) 

By inserting eqs . l |25|26t into the evolution equations (|9) one finds the following recursion relations by combinatorics l Goroff et al. 19861 : 
Jain & Bertschingerlll994l) : 



F<»\k u ...,k n ) = Y - 7 \(2n+ l)a(q 1 ,q 2 )F ( "- m) (k m+1 ,...,k„) + 2 y 6(q 1 ,q 2 )G ( "- m) (k m+1 ,...,k„)l (28) 

*—f (2n + 1) (n - 1) L 1 



m=l 

G^(k u ...,k n ) = y ^ < " l>(k ' , r' k, " ) , [3g(q l ,q 2 )F ( "-"' ) (k m+l ,...,k„) + 2«G ( "-'" ) (k,„ +1 ,...,k„)l , (29) 
(2n + 1) in - 1) L 1 

where qj = k! + ... + k,„, q 2 = k,„ +1 + ... + k„ and F (1) = G (1) = 1. The explicit symmetrised expressions for the second order perturbation 
theory kernels take a very simple and intuitive form: 

- H^HB 2 ^ 

.» (k „ k!) = . (31) 



7 ' 2 hk 2 \k 2 ki) 7 k]k\ 

One can see that mode-coupling to second order reaches its maximum when the contributing modes ki and k 2 are aligned, whereas the kernel 
vanishes for anti-parallel modes. When in eqn. ((26) n different modes qi ...q„ contribute to a mode k, momentum conservation holds, enforced 
by the ^o-distribution: k = qi + ... + q„. Similar to the renormalisation vertex in eqn. (111! also this motivates a diagrammatic description for 
the mode coupling processes (see Fig.QJ. 



3.2 The Four-Point Correlator in Perturbation Theory 

For an analytic expression of the perturbation theory four-point correlator one has to expand the fields in the correlator. Due to the Gaussianity 
of the initial fields ip il) the correlators with an even number of fields ip w will later simplify to products of initial power spectra P l while all 
uneven contributions vanish: 

fa V2 ^ V4 > = <(^ + e" vf + •••) (<P2 } + e" V? + •••) ivf + e" <fif + •••) ivf + & V? + •••)> • (32) 

Simple truncation of the expansion in eqn. M5\ would lead to an inconsistent inclusion of powers of the linear power spectrum P*. We take 
into account all terms up to third order in the linear power spectrum which is equivalent to including terms with initial fields up to sixth order. 
Since the disconnected part will be represented by the full power spectra in eqn. dl4K we are interested in the connected part of the correlator 
only. The connected part can be split into two contributions 

fa V 2 W V4>c = {{ff <P ( 2 2) ) Qpf + all pairs 6 { 1, 2, 3, 4) + (<p™ <pf <pf ^\ + c.p. ( 1, 2, 3, 4} . (33) 

The first contribution originates from second order perturbation theory. In this case two of the fields in the correlator have been expanded to 
second order. The expressions in terms of the initial power spectra and the second order kernels are of the type 

0*1 . k 2 ), (k 3> kO) = 4 e 2 " <5 D (k L .. 4 ) /f/ 5 * 4 ( F®(k 13 , -k 3 ) F< 2) (k 24 , -k 4 ) P[' 3 + F®(k 14 , -k 4 ) F®(k 23 , -k 3 ) P* 14 ) . (34) 
© 2012 RAS, MNRAS OOO.mflol 



6 Gero JUrgens and Matthias Bartelmann 




Figure 1. Diagrammatic representations of the linear propagator g a b(i]a, lb), the linear power spectrum to different times P^(r] a , r\b), the renormalisation vertex 
yit, an d tne perturbation theory kernels F^\q, p) and F^\k, q, p). 



The second contribution is due to third order perturbation theory. Here, one field is expanded to third order while the other three remain at 
linear order. For this reason only one perturbation kernel appears in the expression for this type of contributions: 

t<f (k, , k 2 , k 3 , k 4 ) = 6 e 2 " S D (k L . A ) F®(ki , k 2 , fcj) P*> P* 2 P* 3 . (35) 

With these two functions the connected perturbation theory four-point correlator up to third order in the linear power spectrum P^ can be 
expressed by the following two tree-level contributions 

(ipi V 2 w V 4>c = k2) ' (k3 ' k4) } + a11 P airs 6 !1 ' 2 ' 3 ' 41 + k 2' k 3' k 4> + cp.{l, 2, 3, 4} . (36) 



3.3 Trispectrum in TRG 

Our main objective is to investigate the influence of the perturbation theory trispectrum on the evolution of the power spectrum. Writing the 
connected four-point correlator in terms of the trispectrum T^^l, 

(tpi tp 2 <P3 ^4>c = <5 D (kl... 4 ) TalwW ' ( 37 ) 

we can now include the corresponding corrections into our formalism. Taking the trispectrum in the hierarchy of eqn. (1 1 3 I t into account will 
change the closed system of eqn. d!6t to 

V* = - n ac P* - n fe P k ac + e" J d\ [y%*B%* + (a^ b)] 



■•-"I t iff + y *f 7\ p : k : r + 7 f f " pl r^Y . (38) 



Since the trispectrum is taken from perturbation theory, the evolution of the next higher correlator is not needed for its description. One can 
stick to the same numerical solving procedure as it was presented in section {23}. Modifications appear in the time evolution of the integrals 
^Ldbef f rom e 1 n - G3D v ' a changing the correction integrals A^ cdbe ^ from eqn. <20t : 

K,cd,bef ~* ^-acdjxf + ^Ld,bef (39) 

^U f = ^J^qf^r{^[y^T^ + y^T^ + y^T^ + ( q ^p)]}. (40) 

While the former expression in eqn. ( 120) was a one loop integral only, we now have to integrate twice over the full £-space. The reason for 
this are the additional (So-functions in the disonnected parts of the four-point correlator in eqn. d!4t . The int egration is p erformed numerically 
using Monte Carlo integration techniques from the multidimensional numerical integration library CUBA (Hahn 2005). 
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b a 




Figure 2. Two loop corrections to the power spectrum originating from the trispectrum of tree-level perturbation theory in eqn. ( 1361 . The first three diagrams 
are due to second order perturbations from eqn. (1341 . Third order perturbation theory in eqn. ( 1351 leads to the last diagram. While the perturbation theory was 
included on tree-level only, the non-perturbative time renormalisation leads to 2-loop corrections to the power spectrum. All diagrams are third order in the 
initial power spectrum and second order in the renormalisation vertex 0(y). The latter is due to the fact that both evolution equations in U3t are first order in 
the vertex. 



4 DIAGRAMMATIC DESCRIPTION 

An analytic solution for the system in eqn. d!6t can be formulated ( Pietronil 2008 ). This is still the case for the system in eqn. d38b with 
additional trispectrum terms. Solving first the lin earised evolution equations 1 9) one can write down the linea r solutions for the fields tp at i_(k, rf) 
with the help of the linear propagator g a i,(ri, rf) iMatarrese & Pietronill2007 ; Crocce & Scocci marro 2006b ) : 

Va,h(k, rf) = gab(r], i) fbx(K i) ■ (41) 
Furthermore, the linear propagator has the following properties: 

gabtn, T}') = Sal, 

gab(n,v') gbc(n',n") = &>c(W)- (42) 

With help of this linear propagator a formal analytic solution can be given for the system ((38): 



M = gaAn, 0) g„M o) P k Jn = 0) + J" drftf J d 3 q gjj,, i)g b M i) [y^B^M) + r 



(43) 



gaM 0) g„e(V, 0) g cf (ri, 0) B%f<n = 0) 



• JT d,'e"' fd 3 q[ 2 (y%f P« P'^Tf) + y$> Pfcrf) P*tf) + P^rf) P^rf)) 



r p ' k ' r -i- vPffc-d 

' fdgh Yfgh 



i k.q.r 



f) 



(44) 



For a better understanding of the tripectrum corrections to the power spectrum it is useful to analyse the equations in a diagrammatic 
representation. In the left panel of Fig. Q] symbols for the linear propagator g a b(^la, ??*), the renormalization vertex y^ kb * t and the linear 
power spectrum to different times P\_{r] a , rjb) are depicted. A diagrammatic representation for the perturbation theory kernels are also needed 
to describe the trispectrum. Since we used the perturbation theory trispectrum to third order in the power spectrum, the kernel of second and 
third order are sufficient (Fig.[T}. As one can see from eqs. ( lilt and d26t momentum conservation holds at each intersection in the diagrams. 

The perturbation theory trispectrum terms lead to corrections in the bispectrum in eqn. ( 144b . In eqn. d43t these terms are transported 
to corrections to the power spectrum itself. The types of corrections to the power spectrum are shown in Fig. |2] All the terms are second 
order corrections in the renormalisation vertex 0(y 2 ) since both evolution equations are first order in y. The first three diagrams in Fig. [2] 
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N-body simulations 
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Figure 3. Matter power spectra P* divided by a lineal' spectrum without baryonic wi ggles Pf (Eisenstein & Hull99^) in the mildly non-linear regime. The 
linear spetram is depicted as solid (black) line. A spectrum from JV-body simulations I Carlson et al. 2009) serves as reference (yellow). The dashed (blue) line 
is time-RG theory and the dashed (green) line is time-RG theory with trispectrum included. Left panel: z = 1, right panel: z = 0. 



originate from the second order kernels in eqn. d34| (, while the last diagramm is due to the third order perturbation kernel in eqn. {35} . 
Also here, one can see that all perturbative corrections we included are of third order in the linear power spectrum. While the perturbation 
theory trispectrum is calculated at tree-level, both the evolution equations - for the bispectrum and for the power spectrum - introduce one 
momentum integration. Therefore, all these correction diagrams are two loop diagrams. 

We want to emphazise at this point that the inclusion of the perturbation theory trispectrum does not lead to a simple perturbative 
correction only. At each time step the perturbative trispectrum corrects the evolution of bispectrum and power spectrum. Therefore, from 
that moment on these corrections will be involved in the non-perturbative method of time renormalisation. In this work we only discuss the 
trispectrum correc tions to this me thod, since the quality and performance of the original time renormalisation technique has been thoroughly 
discussed already JPietronill2008l) . 



5 NUMERICAL RESULTS 

We solved the system from equation d38t starting from redshift z = 100 well within the linear regime with linear initial power spectra (see 
eqn. d21t ) and vanishing bispectrum. We evolved the system to redshift z = 1 and z = w ith and without tris pectrum included and compared 
the power spectra with results from numerical simulations of the same ACDM cosmology dCarlson et al.l2009l) . Due to numerical complexity 
only power spectra up to k = 1 /iMpc~' were included in the trispectrum integrals A/j^ . from eqn. ( 140b . However, in the integrals I^dbef 
from eqn. d20^ modes up to k = 10/iMpc" 1 were taken into account, were the results saturate to percent accuracy. 

The results are shown in Fig. [3] in which also the linearly evolved po wer spectrum is depict ed. All spectra were divided by a linear 
power spectrum without wiggles from baryonic acoustic oscillations dEisenstein & Hu 19981 eqn. 29). For z = 1 the results with the 



trispectrum included are in excellent agreement with numerical simulations up to k ^ 0.17/jMpc" 1 . For 0.17/?Mpc~' < k < 0.25 /jMpc~' 
the method performs still better than without trispectrum included. Beyond this regime the perturbative description of the trispectrum breaks 
down and the results are in strong disagreement with simulations. 

It is not surprising that below a certain scale the method performs better without the perturbative trispectrum included. For larger k in 
the integrals AA* crfie ^. in eqn. \40\ also trispectra at smaller scales contribute to the evolution. Since in perturbation theory loop corrections 
become more and more important at smaller scales, the tree-level trispectrum description breaks down above a specific wave number. For this 
reason, beyond k 0.25 /?Mpc~' time renormalisation without trispectrum will lead to better results in comparisons to numerical simulations. 

At z = trispectrum corrections overcompensate the too strong growth in the pure TRG approach on scales below k » 0.15 /?Mpc _1 and 
lead to too little growth. In this regime our results agree with the numerical simulations within 2.5%. The better agreement for z = 1 on these 
scales may simply be due to the breakdown of the tree-level perturbative descrition of the trispectrum at later times. Beyond k » 0.15 /7Mpc~' 
the inclusion of the trispectrum leads to a better agreement with simulations than pure TRG, while both methods show too strong growth 
further inside the non-linear regime. 



The results for pure TRG may differ from the results of Pietroni (2008), since only 12 instead of 14 independent integrals / 



included in the original method. However, in later applications all 14 independent were taken into account. 



acd,bef 
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6 SUMMARY 



In this work the influence of the tree-level trispectrum on the time renormalisation method (Pietroni 2008) has been studied. To keep the 
computational effort at a reasonable level we did not include the full trispectrum evolution from the hierarchy. Instead, we used tree-level 
perturbation theory for the trispectrum in the evolution of the bispectrum. 

(1) The trispetrum was formulated in standard perturbation theory to third order in the linear power spectrum P*. This was included in 
the evolution equation of the bispectrum. 

(2) The numerical method for solving the system was extended by the trispectrum corrections. Due to the linear time dependence of the 
perturbation theory trispectrum it is sufficient to calculate the correction integals at one fixed initial time. Once the corrections are derived 
the extended method operates at the same speed as the time renormalisation method without trispectrum included. The additional 2-loop 
corrections to the power spectrum are of second order in the renormalisation vertex 0(y 2 ) and third order in the linear power spectrum O(P^). 

(3) Perturbative trispectrum corrections are fed into the systems at all times. Once a correction has been included its evolution will 
be described by the non-perturbative formalism of time renormalisation. Therefore, although the trispectrum was only taken into account 
perturbatively, its inclusion can not be interpreted as a pure perturbative correction disentangled from renormalisation. 

(4) We solved the system numerically starting with Gaussian i nitial conditions and a linear power spectrum at an initial redshift of 
z = 100. In comparison to numerical simulations dCarlson e t al. 2009) the inclusions of the trispectrum generally improves the results up to 
k x 0.25 /zMpc~'. However, on larger scales the damping due to perturbative tree-level trispectrum overcompensates the deviation of pure 
TRG from N-body simulations. The results agree with the simulations within 1% up to k as 0.18 /?Mpc~' for z = 1 and within 2.5% up to 
k*0.2hMpc- [ forz = 0. 

(5) Beyond k as 0.25 /iMpc~' the perturbative description of the trispectrum breaks down and the method performs better without 
trispectrum included. This is due to the fact that loop corrections to the trispectrum are not included in our method and become more and 
more important on smaller scales. Adding the perturbation theory trispectrum therefore predominantly pays off on large scales and at the 
beginning of the mildly non-linear regime. 

(6) Although the prediction of the amplitude and position of the first two peaks in the baryonic acoustic oscillations was improved by 
the trispectrum, we are far from reaching 1% accuracy over the entire BAO regime. Also the speed of this method was reduced by including 
the trispectrum, since at a specific time the 2-loop trispectrum corrections have to be derived. Finding a better analytical estimator for the 
trispectrum from other renormalisation approaches, which includes higher order corrections, could improve the results of this method further 
inside the mildly non-linear regime. 
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